First-Principles Study of Electronic and Vibrational Properties of BaHfN 2 
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The transition metal nitride BaHfN2, which consists of weakly bonded neutral slabs of closed 
shell ions, has structural and chemical similarities to other layered nitrides which have impressive 
superconducting T c when electron doped: A^HfNCl, A^ZrNCl, A-TiNCl, with T c - 25.5, 15.2 
and 16.5 K, respectively for appropriate donor (A) concentrations x. These similarities suggest 
the possibility of BaHfN2 being another relatively high T c nitride upon doping, with effects of 
structure and the role of specific transition metal ions yet to be understood. We report first- 
principles electronic structure calculations for stoichiometric BaHfN2 using density functional theory 
with plane-wave basis sets and separable dual-space Gaussian pseudopotentials. An indirect band 
gap of 0.8 eV was obtained and the lowest conduction band is primarily of Hf 5d xy character, similar 
to /3-ZrNCl and a-TiNCl. The two N sites, one in the Hf layer and another one in the Ba layer, were 
found to have very anisotropic Born effective charges (BEC) deviations from the formal charge (-3) 
are opposite for the two sites, and opposite for the two orientations (in-plane, out of plane). LO-TO 
splittings and comparison of BECs and dielectric constant tensors to those of related compounds 
are discussed, and the effect of electron doping on the zone-center phonons is reported. 
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I. INTRODUCTION 

High temperature superconductivity has been a 
puzzle since the quasi-two-dimensional, doped insu- 
lating copper oxides were reported to become su- 
perconducting with very high T c 's. Since then sev- 
eral other layered transition metal oxides have been 
found to be good superconductors although at rela- 
tively low temperature, for example Li^NbO^ and 
Na^CoO^ at about 5 K. The undoped parent com- 
pounds of the cuprate high-temperature supercon- 
ductors are magnetic insulators and their transition 
from a magnetic insulator to a metal upon dop- 
ing completely modifies their electronic structure^. 
These transition metal oxides still attract a great 
deal of interest because the superconductivity is not 
yet well understood. 

Recently, interest has been growing for another 
class of layered superconductors, the transition 
metal nitrides^ such as MNX (M=Ti, Zr, Hf; 
X=C1, Br, I) and ternary transition metal dinitrides 
AMN 2 (A=alkaline earth metal, M=Ti, Zr, Hf), 
some of which have been reported to become su- 
perconducting with high T c values. Superconduc- 
tivity up to 12 K was first measured in /3-ZrNCl 
by Yamanaka et al. in 1996^, and since then the 
highest T c 's that have been measured for these 
transition metal nitrides are as follows: 25.5 K 
for Lio.4 8 (THF) !/ HfNCli, 15.2 K for intercalated /?- 
ZrNC£ and 16.5 K recently reported! for a-TiNCl 
upon doping with Li. These electron-doped transi- 
tion metal nitrides form a new and seemingly uncon- 
ventional class of high T c superconductors because, 
unlike the transition metal oxides, the parent com- 
pounds are not Mott insulators. The parent com- 



pounds for these layered quasi-2D nitrides are non- 
magnetic ionic band insulators with a gap in the 
range of 2-4 eV4"i£ 

In these transition metal nitrides, the supercon- 
ducting mechanism presents a real conundrum. Ex- 
perimental measurements of the isotope effect^ in 
LLjZrNCl show a very weak dependence on the 
N mass, suggesting that electron-phonon mediated 
pairing cannot adequately account for the supercon- 
ductivity in the MNX family. Specific heat measure- 
ments on LLjZrNCli^ estimate the upper limit for 
the electron-phonon coupling constant A ps 0.2 for 
Li .i 2 ZrNCl. A theoretical study on Li x ZrNCli£ pre- 
dicted the coupling constant on the average around 
0.5. The computed and the estimated coupling con- 
stant is far too small to account for the T c of 12-15 
K. The magnetic susceptibility measurements" also 
give a mass enhancement factor that appears too 
small for electron-phonon coupling. 

There is no clear evidence of strong electronic cor- 
relations in the transition metal nitrides. These 
compounds do not show the antiferromagnetism 
that is characteristic of strong correlations nor even 
the Curie- Weiss susceptibility that signals local mo- 
ments, and there is no frustration on cither the hon- 
eycomb lattice or rectangular lattice. In fact, the 
bandwidths of the d states (where doped electrons 
reside) are rather largei^ and undoped systems are 
in d° configurations, so these systems should be well 
described (except for the value of the gap) by first- 
principle calculations employing the local density 
approximation (LDA). There is no observation of 
magnetism in the parent compounds at all, so the 
possibility of spin fluctuations as a pairing mecha- 
nism, similar to what is thought by some to cause su- 



perconductivity in cuprates, seems unlikely. While 
a possible pairing mechanism mediated by magnetic 
fluctuations has been suggested^ this mechanism 
seems at odds with observed behavior so far. Some 
groups have also proposed charge fluctuation o 17 ' 18 
as a pairing mechanism or plasmon enhancement of 
weak BCS superconductivity. 

In this paper we focus on the ternary nitride 
BaHfN2, whose electronic structure and vibrational 
properties have not yet been studied theoretically. 
This compound has many chemical and structural 
similarities with the layered transition metal nitrides 
MNCl's(M=Ti,Hf,Zr) that are impressive supercon- 
ductors when they are electron doped. We suggest 
that this compound has the potential to provide an- 
other high T c transition metal nitride superconduc- 
tor when electron doped. In this paper we present 
calculations of the electronic structure, lattice vibra- 
tions and dielectric constant tensors of BaHfN2, and 
compare them with those of other nitrides. This 
comparison with other layered nitrides may help 
in predicting the origin of the superconductivity in 
these layered nitrides. 

The rest of the paper is organized as follows. 
We first describe the crystal structure of BaHfN 2 
(Sec. |nj and the computational methods (Sec lIII[) . 
Then we present our results for structural and elec- 
tronic properties of BaHfN 2 (Sec. I IV]) . followed by 
analysis of the vibrational properties (Sec. [V]). Fi- 
nally, we discuss the case when BaHfN2 is electron 
doped by replacing one of the two Ba atoms in the 
unit cell by a La atom (Sec. IVI[) . A summary of our 
findings in Sec. IVIII concludes the paper. 



II. STRUCTURE 

We use structural coordinates for BaHfN2 from 
powder X-ray diffraction measurements reported by 
Gregory et al£- The nitrodohafnate BaHfN2 crys- 
tallizes in the tetragonal space group P4/nmm, for 
which KC0O2 is the customary example. The mea- 
sured lattice constants are a = 4.128 A and c = 
8.382 A. There are two inequivalcnt N sites which 
we denote Nl (lying nearly in the Hf plane) and N2 
(nearly in the Ba plane), i.e. BaN2-HfNl. Ba, Hf 
and N2 occupy Wyckoff position 2c (|, i, z) and Nl 
occupies position 2b (|, j, 5) (see Table I). 

The structure of BaHfN 2 shown in Fig. [T]is some- 
times described as composed of [HfN2] 2 ~ anions and 
Ba 2+ cations^. These anions are composed of Hf 
atoms inside a square pyramid of 5 nitrogen atoms 
(four Nl's and one N2), forming layers of edge shar- 
ing pyramids stacked along the c axis; the apices of 
adjoining pyramids arc aligned alternatively up and 
down along the c axis. The Ba 2+ cations are sit- 



uated between these Hf-N sheets levelled with the 
pyramid apices. 

Another description for the structure of BaHfN2 
comprises of (nearly) coplanar (BaN2)~ layer and a 
corrugated (HfNl) + layer, with Hf ions lying alter- 
nately above and below the N2 layer. Each neutral 
BaN2-HfNl-NlHf-N2Ba structural unit (outlined in 
Fig. H]i) is weakly bonded to neighboring units in 
the c direction; upon intercalation, dopant ions will 
lie between the Ba-N2 layers. Similarly, in a-TiNCl 
(Fig. QJ>) there is a layer of transition metal (Ti) 
atoms and N atoms, with Cl _ ions playing a role 
analogous to the (BaN2)~ unit in BaHfN2. Neutral 
TiNCl slabs are weakly bound to each other. A more 
detailed study of that compound will be presented 
separately 

As we show below, the site energies and bond- 
ing of the Nl and N2 sites are quite different, and 
the simple viewpoint of formal closed shell ions, e.g. 
Ba 2+ , Hf 4+ , N 3 ~ may not be sufficient to address 
the structural relationships. 



III. DESCRIPTION OF CALCULATIONS 

We carried out density functional theory (DFT) 
calculations with the ABINIT package^, within 
both local density approximation (LDA) and gradi- 
ent corrected (GGA/PBE 2 ^) exchange and correla- 
tion functionals. Norm-conserving pseudopotentials 
(PSPs) in the relativistic separable dual-space Gaus- 
sian Hartwigsen-Goedecker-Hutter (HGH) for m 19 ' 20 
were used to treat the electronic configuration of Ba 
(5s, 5p, 6s), Hf (5s, 5p, 5d, 6s) and N (2s, 2p), in- 
cluding 5s and 5p semicore states for Ba and Hf. 
Plane wave basis sets with a kinetic energy cutoff 
of 120 Ry were used. A 8 x 8 x 4 Monkorst-Pack^i 
fc-point grids were used to sample the Brillouin zone 
for ground state calculations. We have checked that 
further increasing the cutoff energy to 140 Ry or 
the /c-point grid to 12 x 12 x 6 and 18 x 18 x 6 
has a negligible influence on the relaxed geometry 
and phonon frequencies. The computed LDA band 
structure of BaHfN2 at experimental geometry was 
found to agree well with results obtained from the 
full-potential, all electron code FPLO 2 ^; the latter 
also provides a convenient way to compute contri- 
butions to the electronic bands and density of states 
from individual atomic orbitals. We have confirmed 
that spin-orbit coupling has no significant influence 
to the band structure. Phonon calculations are car- 
ried out at the T point, and the obtained frequencies 
and displacement eigenmodes were used to compute 
Born effective charges and the static dielectric tensor 
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FIG. 1. (Color online) (a) Layered geometry of BaHfN2: N 3 anions form a square base pyramid around Hf 4+ , 
forming a [Hflsb] 2- anion. Ba 2+ cations sit in between the tips of the pyramids, (b) Layered Q-TiNCl structure. 



TABLE I. Structurally optimized lattice constants (a and c) and reduced internal coordinates (z) for atoms in the 
unit cell, obtained from different sets of pseudoptentials (PSPs) including Hartwigsen-Goedecker-Hutter (HGH) 19 ' 20 
and Troullier-Martin (TM)2i type PSPs. The inclusion of semicore 5s and 5p states for Ba and Hf PSPs are indicated 
by Ba sc and Hf sc respectively. The computed band gap for the experimental geometry, Eg Xp , is compared with FPLO 
(FLAPW) results. 





PSP 


Oif &exp 


Of C eX p 


z(Ba) 


*(Hf) 


z(N2) 


E e g xp (eV) 


LDA 


FPLO 
FLAPW 


0.992 


0.986 


0.846 


0.415 


0.177 


0.68 
0.8 




HGH, Ba sc , Hf sc 


0.986 


0.982 


0.848 


0.415 


0.176 


0.82 




TM 


1.005 


1.017 


0.846 


0.415 


0.177 


0.78 




TM, Ba sc 


0.997 


0.987 


0.849 


0.413 


0.170 


0.87 




TM, Ba sc , Hf sc 


0.984 


0.972 


0.847 


0.415 


0.177 


1.11 


PBE 


HGH, Ba sc , Hf sc 


1.000 


1.007 


0.848 


0.415 


0.180 


0.95 




TM 


1.022 


1.036 


0.845 


0.416 


0.186 


1.00 




TM, Ba sc 


1.011 


1.008 


0.850 


0.413 


0.175 


1.13 




TM, Ba sc , Hf sc 


1.004 


0.999 


0.849 


0.414 


0.179 


1.25 


Exp.' 








0.8479 


0.4142 


0.168 





culated structural, electronic and vibrational prop- 
erties of BaHfN2, especially the inclusion of 5s, 
5p semicore states for Ba and Hf, we also carried 
DFT calculations using Troullier-Martin (TM)2i 
PSPs generated using the fhi98PP program^ with 
LDA/PZ ZL and GGA/PBE exchange-correlation 
functionals, respectively. A 8 x 8 x 4 fc-grid and 
a kinetic energy cut off of 90 Ry were used. 

To simulate the doped BaHfN2, we replaced one 
of the two Ba atoms in the unit cell with La, which 
provides an extra electron per formula unit and pro- 
vides metallic screening with its impact on the zone- 
center phonons. The HGH PSPs were used to do 
calculations for the doped system with kinetic en- 
ergy cut-off of 120 Ry and and a k-point mesh of 
12 x 12 x 6. 



IV. ELECTRONIC STRUCTURE 
CALCULATIONS 

A. Structural Relaxation and Electronic 
Structure 

It is instructive for future studies to quantify the 
effects of different types of pseudopotcntials on the 
relaxed structure and the energy gap. We have per- 
formed structural relaxation for the BaHfN2 crys- 
tal using both LDA and GGA/PBE, and results are 
compared with experimental geometry^ in Table HI 
In addition to the HGH PSPs, we also employed 
different sets of norm-conserving TM PSPs. The in- 
clusion of 5s and 5p semicore states, if present, is 
denoted by Ba sc and Hf sc respectively. The opti- 
mized cell parameters and internal coordinates were 



compared with all-electron, full potential reference 
results from FPLO^ The band gap was computed 
for the experimental geometry, E° xp using FPLO^ 
and full potential linear augmented plane waves code 
(FP-LAPW) implemented in the Elk code^ The re- 
sults are listed in Table |U 

One observation from Table U is that as one in- 
cludes semicore states of Ba and Hf, the equilibrium 
lattice constants becomes smaller, while the energy 
gap Eg Xp increases substantially. Such trend holds 
for both LDA and PBE PSPs of TM type. The in- 
clusion of semi-core states is important since there 
is significant amount of hybridization of Ba semi- 
core 5p states with N2 2s state, as can be seen from 
Fig. [2] We observed a trend of decrease in equilib- 
rium lattice constants upon inclusion of semi-core 
states. This is mainly due to the decrease of Ba-N2 
and Ba-Hf bond length. 

On the other hand, TM-typc PSPs generated us- 
ing the fhi98PP code are so-called single projector 
pscudopotcntials, e.g. there is only one pscudopo- 
tential for each angular momentum type, not for 
each valence orbital. The corresponding PSPs with 
Ba sc and/or Hf sc failed to describe the energy posi- 
tion of 6s states properly, even for the isolated atoms. 
Such PSPs tend to predict energy gaps of BaHfN2 
larger than those without the semicore stats in the 
valence configuration. The discrepancy in E° xp com- 
pared with all-electron calculations can be as much 
as 0.4 eV. The HGH pscudopotcntials, on the other 
hand, were constructed with multiple projectors per 
angular momentum type and therefore can describe 
orbitals of same angular momentum but different 
shells reasonably well. Indeed, from Table U we find 
that overall HGH pscudopotcntials give structural 
properties and energy gap similar to those calculated 
from all-electron calculation and in the following, we 
present results obtained with LDA type HGH pseu- 
dopotentials. 

Figure [3] shows the band structure of BaHfN2, us- 
ing the so-called fatbands emphasis of band charac- 
ter for Hf 5d xy states. The fatbands are obtained by 
using the expansion of the wavefunctions in terms of 
the basis atomic orbitals at each k-point 
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FIG. 2. (Color online) Projected density of 

States(PDOS) showing strong hybridization of Ba 5p 
semi-core states with N2 2s states and a weak mixing 
with N2 2p states. 




FIG. 3. (a) Band structure showing Hf 5d bands. The 
amount of Hf 5d xy character is shown in the band struc- 
ture by the width of the lines, revealing that the lowest 
conduction band has strong Hf 5d xy character, (b) The 
density of states, units for horizontal axis are eV -1 . The 
shaded area indicates the amount of Hf 5d character. 
The Fermi level lies at the bottom of the lowest conduc- 
tion band. 
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where n is the band index and L = Im is the orbital 
index. R + s denotes a regular lattice site, with R 
a Bravais lattice vector and s a basis vector of the 
unit cell. The width of the fatband is proportional 
to |c^| 2 . 

BaHfN2 is a band insulator with a calculated band 
gap 0.68-0.80 cV using full-potential, all-electron 
methods (FPLO, LAPW). Given the usual LDA un- 
derestimate of band gaps the true gap of BaHfN 2 



may be as large as 1.5 eV. This layered ionic semi- 
conductor character is very similar to that of the 
MNC1 compounds (M=Ti, Hf, Zr) which have been 
found to superconduct with impressively high T c val- 
ues when electron-doped. The lowest conduction 
band in BaHfN2 has primarily Hf 5d xy character 
with a width of 3 eV. Since these states are empty, Hf 
is formally 4+, and the rest of the electronic struc- 
ture is indicative of a closed shell, ionic insulator 
with some mixing of N 2p states and Hf hd states. 
This characterization is similar to ZrNCl, which is 




5 

Energy (eV) 



in 



FIG. 4. (Color online) Total density of states (DOS) for 
BaHfN2, with projected density of states for (a) Ba, Hf 
with (b) blowup of conduction bands, and (c) Nl, N2 
with (d) blowup of conduction bands. The lowest con- 
duction bands, are composed of Hf 5d xy states, (b) and 
(d) show almost 2D like behavior around the bottom of 
the conduction band up to ~0.9 eV, where Ba conduc- 
tion states start to appear. 



also an ionic semiconductor with lowest conduction 
band having Zr in-plane Ad character^ and TiNCl 
which has Ti 3d xy character @|29J]. The Hf 5d char- 
acter extends through a range of 8 eV beyond the 
Fermi level (see Fig. 2b) partially due to crystal field 
splitting of the 5d orbitals. 

Figure [4] shows the projected density of states in 
the valence-conduction band region. Integrating the 
density of states, we find that only above a dop- 
ing level of 0.17 electrons (shifting the Fermi level 
upward by 0.9 eV) do bands other than the two- 
dimensional Hf 5d xy band start filling up with elec- 
trons. The conduction bands which appear at 0.9 
eV above the Fermi level are Ba 5d states at M that 
have some dispersion along k z (compare the points 
M and A) . There is also some nitrogen hybridization 
in the conduction bands that contribute to their in- 
plane dispersion. 

As mentioned earlier, there are two N sites, cor- 
responding to the corrugated layers Hf 4+ Nl 3_ and 
the nearly flat Ba 2+ N2 3_ layers. The projected DOS 
shows that the Nl and N2 ions are quite distinct in 
electronic character. The N2 ion in the BaN2 layer 
has the more weakly bound (and therefore more po- 
larizable) N 2p states, lying just below the gap. The 
energy states of Nl ion in the HfNl layer (closer 
to the highly charged Hf 4+ ion) is centered about 2 
eV lower in energy. From purely energetic (bind- 
ing) consideration, Nl should have a correspond- 
ingly lower polarizability. 



B. Fermi Surface for Electron Doping 

According to our band structure calculations in 
FigO when BaHfN2 is electron-doped (as are ZrNCl, 
HfNCl, and TiNCl when they become superconduct- 
ing), the Fermi surface will be a single r-centercd 
(nearly) circular surface up to x=0.17 for doping 
concentration x when doped with alkali metals; for 
doping concentration x > 0.17 carriers also go into 
the bottom of the Ba 5d bands at M and then addi- 
tional Hf 5d bands at T. In this respect, BaHfN2 is 
similar to TiNCl^S, but different from the hexagonal 
compounds [(Zr,Hf)NCl] which have the band mini- 
mum at the zone corner K points, of which there are 
two. 

This difference in Fermi surfaces has some im- 
portance for electronic response. In the doped 
(Zr,Hf)NCl compounds, nesting of the two Fermi 
surfaces has recently been put forward 3 ^ as a poten- 
tial source of spin fluctuations, which was suggested 
as a possible candidate for pairing mechanism. How- 
ever, with a single simple Fermi surface such as dis- 
played by doped TiNCl (known to be an excellent 
superconductor) and doped BaHfN 2 (which we sug- 
gest by analogy may be a good superconductor) , this 
mechanism is not available. Since the superconduct- 
ing T c is large (and almost similar in magnitude) in 
TiNCl and ZrNCl, and their characters are other- 
wise so similar but Fermi surfaces are different, the 
mechanism of spin fluctuations seems to be degraded 
in likelihood. Electron-phonon coupling is weak in 
A x ZrNCl, and the materials are pauli-paramagnctic. 
Thus for possible pairing in these materials, long 
sought electronic mechanisms of pairing need con- 
sideration. A single electron gas with a given value 
of kp is different from a pair of identical, degenerate 
electron gases with a value of kp that is \j\[2 as 
large. One main difference, as mentioned above, is 
that there is no nesting that might enhance charge 
fluctuations (as is the case also for spin fluctuations) 
or affect pairing symmetry. Another clear difference 
is that the characteristic momentum scale kp is dif- 
ferent. 



V. VIBRATIONAL SPECTRUM 

A. Born-Effective Charges 

The MNC1 compounds become superconducting 
upon electron doping from an ionic insulator to the 
metallic phase. The relatively low density electron 
gas that is formed upon light doping might not ad- 
equately screen the ionic nature of the MN layers, 
so the electronic response may still have short-range 
ionic character. For this reason we calculate and 



TABLE II. Born effective charges (BECs) of BaHfN 2 and metallochloronitrides MNC1 (M=Hf,Zr,Ti). All BECs 
have been calculated at experimental lattice constants and relaxed atomic positions. a-TiNCl is orthorhombic, so 



7^ Z* y and both are listed. 



For comparison, BECs of NaCoOz"*- are also listed. 



7* 17* 

^xx I ^yy 

yFormal 



BaHfN, 



Hf Nl N2 Ba 
4.52 -4.66 -2.59 2.73 
3.09 -1.65 -4.58 3.14 

+4 -3 -3 +2 



/3-HfNCl 



Hf N CI 
4.7 -3.4 -1.2 
2.5 -1.6 -0.9 
+4 -3 -1 



/3-ZrNCl 



Zr N CI 
5.1 -3.8 -1.3 
2.6 -1.6 -1.0 
+4 -3 -1 



a-TiNCl 



Ti 

5.9/6.4 

1.7 

+4 



N CI 

■5.5/-4.4 -0.4/-1.9 

-0.9 -0.8 

-3 -1 



NaCo0 2 



Na Co O 
0.87 2.49 -1.68 
1.37 0.87 -1.12 
+1 +3 -2 



analyze the Born effective charges in some detail. 

The Born effective charge (BEC) tensor Z* is a 
fundamental quantity for the study of lattice dynam- 
ics, describing the long range Coulomb part of the 
force constants. The Born effective charge Z* of 
atom k can be viewed either as the change of po- 
larization P 1 induced by the periodic displacement 
r K]Q . under the condition of zero macroscopic elec- 
tric field, or as the force F Ka induced on atom k by 
an electric field £ 7 under the condition of no atomic 
displacement. It can also be expressed as the second 
partial derivative of the total energy with respect to 
the displacement and the electric field: 



K " a dr K .. n 9£v 



d 2 E 






(2) 



where V is the volume of the unit cell. One might 
naively expect the BECs to be close to the formal 
charges of the compound, but this is often not the 
case. The BEC can be decomposed into the charge of 
the (pseudo-) ion k, Z k , and the electronic screening 
term, 



A 7 



7* — 7 



AZ 



K,"fCt 



(3) 



The computed BECs for BaHfN2 are provided in 
Table |TTJ In tetragonal symmetry the BEC tensor 
is diagonal and reduces to two values Z* x = Z* 
and Z* z . The BECs for Hf and Ba in the plane are 
reasonably close to their formal charges (Hf is 0.52 
larger and Ba is 0.73 larger), while their perpen- 
dicular charges differ substantially from the formal 
charges (being smaller than the formal charge for 
Hf) indicating a more complex electronic response. 

The two nitrogen sites have very different 
and unusual BECs. For comparison, previously 
calculated^! BECs of NaCo02 are also listed in Ta- 
ble HU The BECs for O in NaCo0 2 are rather un- 
interesting, both being smaller than their respective 
formal charges. By contrast, in BaHfN2, Nl has a 
BEC of -4.66 in the plane and -1.65 perpendicular to 
the plane, respectively with magnitude much larger 
and much smaller than the formal charge. N2, on the 
other hand, has a BEC of -2.6 in the plane and about 
-4.6 perpendicular to the plane, again very differ- 
ent from the formal charge but in the opposite sense 



with respect to Nl. The BEC of Nl is consistent 
with covalent bonding between Nl and Hf, given its 
anomalously large magnitude in the plane. N2 be- 
haves in the opposite way, and its BEC is consistent 
with little covalent bonding with Ba (as expected) 
but with significant covalent interaction with Hf in 
the inner layer. BECs with magnitudes greater than 
the formal charges reflect large electronic response 
to atomic motion. For example, in the case of pcr- 
ovskitc BaHfOjr^ the large Born-effective charges of 
Hf (Z* = 5.75 ) and O ( Z^ = -4.42, Z* ox = -2.03) 
indicate a mixed ionic-covalent nature of Hf-0 bond, 
similar to the case found here for the Hf-Nl bond. 
Also, Ba in BaHfC>3, which has a cubic site symme- 
try, was found to have a similar average BEC (Z* = 
2.72 as that computed for Ba in BaHfN2 here ( Z* 
= 2.87)). 

The Hf-Nl layer, taken as a unit, behaves as if 
having a charge of -0.14 (nearly neutral) in the plane 
and +1.44 (slightly cationic, rather consistent with 
the formal charges) for vibrations perpendicular to 
the plane. The Ba-N2 layer behaves in an opposite 
manner. 

In Table HH we draw a comparison be- 
tween the BECs of metallochloronitrides MNC1 
(M=Ti,Hf,Zr)22 and those of MNC1 (M=Ti, Hf, 
Zr)22 show similar trends as do the BECs of Nl 
in BaHfN 2 and the N in MNC1. There is consid- 
erable anisotropy in the effective charges for both 
the M and N ions. The effective charge for CI in 
MNC1 is close to its formal charge (reflecting its high 
electroncgitivity), however the CI analog in BaHfN2, 
Ba-N2, is somewhat different, with large anisotropy. 



B. Zone-Center Phonons 

BaHfN2 has 8 atoms in the unit cell resulting in 
24 phonon modes, three of which arc acoustic modes 
and the remaining 21 are optical modes. The phonon 
frequencies are listed in Table IIIII with their polar- 
ization and symmetry. There are 8E g + 6E U modes 
with polarization perpendicular to the c axis (within 



3A 2 



the x-y plane) and 3Ai g 
polarization along the z axis 



lBi a modes with 



'is 



The modes Ai g , B 



Fig are Raman active. 



and the modes A^ Ul E„ 



Iff! 

are 



TABLE III. Calculated zone-center phonon frequencies for BaHflSh and Bao.sLao.sHfN^. The phonons were computed 
using the optimized geometry. All phonons with x-y polarization are doubly-degenerate. The symmetry column refers 
to the symmetry of the phonons in higher symmetry insulating system (point group D^h vs. C± v for the metallic 
system). There is a split in degeneracy in the long wavelength limit between the LO and TO modes, and the 



magnitude of LO-TO splitting, s/L 



LO 



is listed in the last column. All frequencies reported in cm 



Mode Symmetry Polarization BaHfN2 frequency Bao.sLao.sHff^ 









LUTO 


ULO 


/ 2 2 


frequency 


l-2 a 


En 


x-y 


72 


93 


59 


76 


3-4 b 


E fl 


x-y 


82 






94 


5 a 


A 2u 


z 


105 


144 


98 


144 


6 b 


Ai s 


z 


120 






136 


7-8 b 


E fl 


x-y 


152 






148 


9 b 


Ai 9 


z 


172 






175 


10-ll a 


E u 


x-y 


210 


240 


116 


210 


12-13 b 


E fl 


x-y 


232 






283 


14 b 


Big 


z 


341 






328 


15-16 a 


E„ 


x-y 


424 


614 


444 


475 


17 a 


A 2u 


z 


468 


492 


152 


457 


18-19 b 


E, 


x-y 


623 






651 


20 a 


A 2u 


z 


641 


751 


391 


596 


21 b 


Ai 9 


z 


717 






646 



a IR active mode 

b Raman active mode 



infrared active. We found some phonon frequen- 
cies were quite sensitive to the inclusion of semi-core 
states in Ba and Hf pseudopotentials with several 
phonons differing by 15-30% in the absence of semi- 
core states. 

The LO-TO splittings of the IR active modes can 
be related to the Born-effective charges as 



V" 



[2 2 1 

l^LO.m ~ u TO,m\ 



4-7T 

e°° V 



E 



(ez; 



M K 



(4) 



In this relation, m goes over the IR active modes 
of a given polarization direction a, M K is the ionic 
mass of the atom k, and e^ a is the a-th diagonal 
element of the high frequency dielectric constant. 
When the LO phonons are excited, a macroscopic 
electric field is created due to the long range na- 
ture of the Coulomb interaction. The squares of the 
BECs, divided by the mass, give the contribution of 



TABLE IV. Calculated macroscopic dielectric constants 
for BaHfN2 and group IVB nitrochlorides using (a) fully 
relaxed geometry and (b) the experimental structure 
with relaxed atomic positions. 



*lvv 



x/yy 



(a) BaHfN 2 7.47 7.55 33.8 21.4 

(b) BaHfN 2 7.35 7.31 44.7 24.0 
(b) Q-TiNCl 6.9/7.4 3.2 22.3/38.3 3.7 
(b) /3-ZrNCl 6.2 4.4 13.8 5.9 
(b) /3-HfNCl 5.4 4.0 11.1 5.1 



that ion to the electric field. One interesting exam- 
ple (see Table ITO]) : the large splitting (Aw) of mode 
15 is due largely to the fact that the light Nl ions 
are vibrating in the x-y plane and the BEC for Nl is 
rather large (-4.66), accounting for most of the shift 
of 45% in LO frequency. The second largest split- 
ting (Aw) is for the mode 20, with a shift of about 
17% in frequency. In this mode we have primarily 
N2 vibrating along the z direction and the BEC for 
N2 along z is Z* Z (N2) = -4.58 which again accounts 
for the large splitting. Mode 5 has a large relative 
shift of 37%; it has primarily Nl and Ba vibrations 
opposite to each other. 



TABLE V. Lattice contribution (defined in Eq. [6]) to 
the macroscopic dielectric constants for selected phonons 
computed using fully relaxed geometry. Only 3 phonon 
modes contribute to the lattice polarizability in each di- 
rection. 



Mode 


euro (cm i ) 


47r_f77i,a;:E 


47rP m , 22 


1-2 


72 


14.0 




5 


105 




10.1 


10-11 


210 


5.1 




15-16 


424 


7.3 




17 


468 




1.3 


20 


641 




2.4 



C. Dielectric Response 

We now discuss the electronic and ionic lat- 
tice contributions to the macroscopic dielectric con- 
stants, computed for BaHfN2 by doing the phonon 
calculation using both experimental lattice con- 
stants with relaxed geometry and completely relaxed 
geometry. e°° denotes the high frequency electronic 
response where there is no contribution from the 
ionic lattice polarizability (p Mn ) and e° is the sum of 
the electronic and ionic response. Their relationship 
is given by 

e% = e d + AttP™ = e^ + 4^ £ P mM , (5) 

771 

„ 1 Jrn,o>j3 



m,aj3 



n uj 2 ' 



(6) 



where the sum is over all the modes and S mta p is 
the mode-oscillator strength tensor which is defined 
as 

S m ,af3 = ^2 Z taa' U m(K,a')x ^ Z*, >t3f3 ,U m («', /?'), 

(7) 
where U m {K,a) is the component of the phonon 
eigenvector for the m-th mode corresponding to the 
displacement of atom k in direction a. 

The values of the static dielectric constants 
for BaHfN2 and comparison to the group IVB 
metallochloronitrides^ are given in Table IIVI e°° 
for BaHfN2 is larger than for the mctallochloroni- 
trides (the electronic polarizability e°°-l is 50-60% 
larger than the Hf counterpart), consistent with the 
smaller band gap (0.8 eV versus around 1.8 eV for 
ZrNCl and HfNCl). 

With the exception of the in-plane values for 
TiNCl, the lattice polarizability e° — e°° of other 
MNCl's is smaller by a factor of 5-15 relative to that 
of BaHfN2. This is due to the fact that the phonons 
in BaHfN2 are softer than the phonons in MNCl's 
and the oscillator strengths for modes which con- 
tribute to P lon in BaHfN2 are much larger than the 
modes contributing P lon for the other chloronitrides. 
The modes that contribute to lattice polarizability 
P ion in BaHfN 2 are shown in Table ED together with 
their contributions. The main contribution arises 
from the lowest IR active mode for each direction (x 
and z). We found that the frequency of the lowest 
IR mode in the plane is highly sensitive to the inclu- 
sion of the semi-core states in Ba pseudopotential, 
and as a result when these states are absent e xx is 
reduced from 44.7 to 17.3. 

For some comparison we note the dielectric con- 
stants for a few transition metal nitrides. The high 
frequency dielectric constants for group IVB nitrides 
have been reported^ as follows : TisN4, Zr3N4, 



Hf 3 N 4 with e 00 = 18.31, 9.36, 10.10 respectively. 
TisN4 has a higher dielectric constant due to a small 
band gap. 



VI. DOPING WITH ELECTRONS 

We consider La substitutional doping by replac- 
ing one Ba in the unit cell with La. The vibra- 
tional frequencies for the doped system are included 
in Tabic IIII1 Since there is no experimental data 
on Bao.5Lao.5HfN2 we use relaxed lattice constants, 
which arc smaller by almost 3% than the experimen- 
tal lattice constants of BaHfN2. The band struc- 
ture of (BaLa)o.5HfN2 is shown in Fig. [5] with fat 
bands for Hf and La. The nearly dispersionless 
La 4f states are located 1 eV above Ep but have 
no clear impact on what we discuss in the follow- 
ing. There are two nearly cylindrical Fermi surfaces 
around T, one of which has Hf hd xy character and 
another has La 5d x 2_ y 2 character; mixing may oc- 
cur at or near crossing of the Fermi surfaces. The 
lowest conduction bands near M, primarily La 5d yz , 
5d xz in character (not shown as fatbands) arc signif- 
icantly lowered from the corresponding undoped Ba 
bands, creating two Fermi surfaces that are larger 
than anticipated from a rigid band picture using the 
BaHfN2 bands. These surfaces have significant three 
dimensional character, but even at this (large) dop- 
ing level they do not reach the top of the zone (the A 
point). The lowering of the bands around M is also 
seen within the virtual crystal approximation (where 
both Ba and La are replaced by an 'average ion') so 
we expect this to be a robust feature for doping by 
La. A small Fermi surface near M, of mainly La 5d z 2 
character, arises near this level of doping. 

The phonon frequencies for the metallic system 
are shown in Table UITl Although the replacement of 
one Ba with La changes the symmetry, correspond- 
ing modes between the two systems can be identified 
by examining the scalar products of their eigenvec- 
tors. Several of the softer modes have slightly higher 
frequencies in (BaLa)o.5HfN2 likely due to the de- 
creased lattice constant of the metallic system. More 
interestingly, several of the high frequency modes are 
rcnormalized in the metallic system. These high fre- 
quency modes are dominated by motion of the Nl 
and N2 atoms. 



VII. SUMMARY 

We have examined the electronic and vibrational 
structure of the ternary nitride BaHfN2 within den- 
sity functional theory. We find that BaHfN2 seems 
to have chemical and electronic similarities with 




FIG. 5. Doped band structure showing La and Hf 5d 
bands. The amount of Hf 5d xy and La 5d x 2_ y 2 character 
is shown in the band structure by the weight of the points 
in the lines. The lowering of the 5d character due to 
replacement of Ba by La is substantial. The Fermi level 
lies at an energy of eV. 



high T c mctallochloronitrides MNCl's (M=Ti, Hf, 
Zr), so its candidacy as another high T c supercon- 
ducting nitride is plausible. The basic electronic 
and vibrational properties of the undoped insulating 
phase provide a basis for an understanding of the 
behavior of BaHfN2 upon doping. We find highly 
anisotropic Born effective charges for the N ions, 
with anisotropics that have an opposite sign for the 
two N sites. These differences suggest unusual elec- 
tronic screening in BaHfN2 , so that in the absence of 
strong correlations and magnetism, electron-electron 
interactions might play a significant role in pairing 
in these layered nitrides when electron doped. The 
large BECs result in large LO-TO splittings for some 
zone center phonons, as well as large dielectric con- 
stants that also imply unusual characteristics of elec- 
tronic screening. We also provided initial analysis of 
how the system is affected by doping and found that 
bands near M are significantly lowered so that con- 
duction would occur both in the Hf states near T 
and in states near M located in the BaN layer. 

The MNC1 compounds, which are impressive su- 
perconductors when doped, provide an interesting 
analogy to BaHfN2 , with its similar structural, elec- 
tronic, and vibrational similarities but larger elec- 
tronic screening. One potentially important differ- 
ence is worth noting. From the point of view of vapor 
phase growth of the materials, MNC1 contain two 
reactive, highly electronegative anions. There is rel- 
atively little experience in vapor growth (molecular 
beam epitaxy, pulsed laser deposition) of such mate- 
rials. Since the appearance of the high T c cuprates, 
there has been a huge amount of experience accumu- 
lated, and expertise gained, in deposition of oxides 



with one electronegative anion but several cations. 
From this viewpoint, BaHfN2 - with one anion and 
two cations - is attractive for vapor phase growth, 
and hence for study and potential application of 
ultra-thin superconducting layers. 
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Appendix A: Description of the vibrational 
modes 



Here we provide a brief characterization of the 
eigenvectors of all optical modes at q — 0, which can 
be useful in interpreting optical data and in compar- 
ing with similar ionic semiconductors. The units are 
cm -1 . 

uj = 72 : Ba oscillating against other atoms in 
the x-y plane. These two degenerate soft mode is 
largely responsible for the large static dielectric con- 
stant e° / . 

xx/yy 

uj = 82 : This mode is primarily out-of-phase Ba 
vibrations in the x-y plane. 

uj = 105 : Ba and Nl moving against each other 
with strong amplitude along the z axis. Hf and N2 
are in phase with each other oscillating weakly as 
compared to Ba and Nl. 

uj = 120 : This mode is primarily out-of-phase Ba 
vibrations along the z axis. 

uj = 152 : This mode has primarily N2 and Hf 
vibrations in the x-y plane with N2 vibrating with a 
large amplitude as compared to Hf with Ba and Nl 
participating very weakly. 

uj = 172 : Hf and N2 vibrating along z, with Ba 
vibrating very weakly opposite to Hf. 

uj = 210 : These two mode are primarily in-phase 
N2 vibrations in the x-y plane. 

uj = 232 : These two modes are primarily out-of- 
phase N2 vibrations in the x-y plane. 

uj = 341 : This mode exhibits pure Nl vibrations 
along the z axis, with nearest neighbor Nl atoms 
out of phase. 

uj = 424 : These two modes have primarily in- 
phase Nl vibrations in the x-y plane with Hf mov- 
ing weakly against Nl and Ba and N2 in the Ba-N2 
moving weakly in phase with Nl. 

u> = 468 : Nl vibrating with large amplitude in 
phase with Ba, N2 and out of phase with Hf along 



the z axis and each of the other atoms participate 
very weakly. 

uj = 623 : These two modes are primarily out-of- 
phase Nl vibrations in the x-y plane. 

u> = 641 : This mode is primarily in-phase N2 
vibrations against other atoms along the z axis. This 



mode dominates the ionic response in the z direction, 
especially due to the large magnitude of BEC for N2: 
Z* Z (N2) = -4.58. 

uj = 717 : This mode is primarily out-of-phase N2 
vibrations along the z axis with very weak partici- 
pation from Ba and Hf. 
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